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QUASISTATIONARY DISTRIBUTIONS FOR ONE-DIMENSIONAL 
DIFFUSIONS WITH KILLING 

DAVID STEINSALTZ AND STEVEN N. EVANS 



Abstract. We extend some results on the convergence of one-dimensional 
diffusions killed at the boundary, conditioned on extended survival, to the 
case of general killing on the interior. We show, under fairly general conditions, 
that a diffusion conditioned on long survival either runs off to infinity almost 
surely, or almost surely converges to a quasistationary distribution given by 
the lowest eigenfunction of the generator. In the absence of internal killing, 
only a sufficiently strong inward drift can keep the process close to the origin, 
to allow convergence in distribution. An alternative, that arises when general 
h^ ' killing is allowed, is that the conditioned process is held near the origin by a 

high rate of killing near oo. We also extend, to the case of general killing, the 
standard result on convergence to a quasistationary distribution of a diffusion 
C^ ' on a compact interval. 



1. Motivation 



Among the puzzles that occupy mathematical demographers, one of the most 

tantalizing is the phenomenon sometimes called "mortality plateaus" . While the 

^S I knowledge that human mortality rates increase with age during the years of matu- 

^^ . rity remounts to the primitive past of statistical science, only in recent years has it 

^D ' become apparent that this accelerating decrepitude slows in extreme old age, and 

r~| ■ may even stop |Vau97j . Analogous flattening of the mortality curves for fruit flies 

[SLOV92 is well established. 

a"* . In their attempts to explain this widespread, and possibly near-universal phe- 

nomenon, mathematical demographers have turned repeatedly to Markov models 

^ . of mortality. These are stochastic processes where "death" is identified with a ran- 

dom stopping time, which arises from a typically unobserved Markov process. Some 
models, such as the "cascading failures" model of H. Le Bras |Le 761 (described in 

C^ ' section EJ, were originally introduced to model the classical exponentially increas- 

ing "Gompertz" mortality curve, and were only later shown jGG91| to converge to 
a constant plateau mortality rate. Others, such as the series-parallel model of L. 
Gavrilov and N. Gavrilova GG91 , and the drifting Brownian motion model of J. 
Weitz and H. Eraser |WF01| . were introduced explicitly with mortality plateaus in 
mind. (The drifting Brownian motion, though, it should be pointed out, was first 
brought into demography by J. Anderson AndOO with slightly different motives. 
The later work of Weitz and Fraser was apparently independent.) 
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2 DAVID STEINSALTZ AND STEVEN N. EVANS 

While these models illustrate theories of the development of senescence, they 
are not really explanations to be set beside the usual demographic explanations for 
mortality plateaus. The authors derive the convergence to constant mortality from 
explicit computations alone, which cannot help to make palpable the driving force 
behind the convergence. Nor do they suggest whether the convergence is a typical 
phenomenon, or whether it depends on the many arbitrary particularities of the 
models. 

As we have explained in |SE04| . mortality plateaus in Markov aging models are 
examples of the generic property of convergence to quasistationary distributions. 
What is more, translating quasistationary distributions into demographic reality 
offers a novel account of the mortality plateaus, distinct from the standard dyad 
of population heterogeneity (mortality rates stop rising because the few survivors 
at extreme ages are an intrinsically healthier subset of the initial population) and 
temporal heterogeneity (the aging process itself slows down in time) . 

One might have expected, in a field as thoroughly ploughed over as the asymp- 
totic behavior of Markov processes, that any problem thrown up by applied science 
would already have been effectively answered in the probability literature. This 
problem seems to be an exception. In particular, one-dimensional diffusions which 
are killed during their motion appear to have been overlooked. The case of a diffu- 
sion with one inaccessible boundary and one regular boundary, at which it is killed 
(or partly killed and partly reflected), with no internal killing, has been elucidated 
by P. Mandl [Man61[ . His methods were further developed, in the more specialized 
case of pure killing at the boundary, by P. Collet, S. Martinez, P. Picco and J. San 
Martin |CMM95llMP'tiM98llMMmllMMn4J . 

Some of the general principles are the same as those for countable discrete-time 
Markov chains, first discussed in |Yag47| (in the context of branching processes), 
and further developed by D. Vere- Jones alone jVJ62| (where the R-theory, later de- 
veloped by R. Tweedie |Twe74| . was first introduced) and with E. Seneta SVJ66] . 
The corresponding theory for continuous-time countable-state processes was intro- 
duced by J. F. C. Kingman [Kin63j . We note that the results of |CMM95] are pre- 
saged in the work of E. van Doom jDooQlj and of M. Kijima and E. Seneta jKS91| 
on birth-death chains conditioned on not having dropped below 0, which itself is 
built upon a birth-death-chain countable version of Mandl's work by J. Cavender 
|Uav78) . The study of quasistationary distributions for birth-death chains was ini- 
tiated by P. Good |(tOo68| . There is a strong analogy between birth-death chains 
killed at and diffusions killed at 0. Some features of the problem, in particular the 
role of the initial distribution, are easier to understand in the birth-death chains, 
and are well described by |Doo91| in particular. A slight generalization is offered 
by |FKMP95] , where the Markov process in the positive integers is allowed to have 
arbitrary jumps, but is still constrained to have killing only at 0, and to take a 
long time to reach from far away. Similarly, quasistationary distributions for a 
particular density-dependent branching process were studied by G. Hognas |Hog97| . 

More general killing has primarily been considered in the finite state-space set- 
ting. General theory for quasistationary distributions on finite state-spaces is 
worked out in "DSGS", and at greater length in |Sen73| : applications to genetics 
are discussed in Sen66 . In this context, it is not hard to show that quasistation- 
ary distributions exist (under basic irreducibility assumptions), and also that the 
distribution of the process conditioned on survival converges to that distribution 
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(sometimes referred to as the "Yaglom limit"). Only recently have comparable 
results been developed for Markov chains on infinite state spaces, by F. Gosselin 
[GosOl" (for countable state spaces) and by J. B. Lasserre and C. E. M. Pearce 
[LPOI (for general state spaces) . Both of these have results similar to those of the 
present work (allowing for the substantial technical difference in setting), but also 
fundamentally different, in that they rely on a Lyapunov-type criterion. We are 
not aware of any earlier work that derives these kinds of convergence results for 
diffusions with internal killing. 

In the present work, we extend the results of jManGlj and JCMM95' to the case 
of general killing, and prove new results which would not arise when killing occurs 
only at 0. The essential result is a dichotomy: The process either escapes to oo with 
probability 1 (so that the probability of finding it in any compact interval goes to 0), 
or it settles into the quasistationary distribution given by the top eigenfunction of 
the adjoint generator. Along the way, we shore up the foundations of this subject, 
supplying details which were left vague. 

Our methods are more "hands-on" than those of JCMM95] . In part, this is 
compelled by a serious gap in their proof. In their Lemma 2 it is stated that 
a sequence t„ — > oo may be found such that the measures ¥fi{Xt^ & ■ \ Tg > tnj 
converge, where h is an initial distribution with compactly supported density. Their 
proof simply attributes this fact to compactness of the space of probability measures 
on the extended real half-line. While this is, of course, correct for any single h, they 
claim, with no further justification, that the same limit measure may be obtained for 
an arbitrary countable set of initial distributions. This is equivalent, in the notation 
that we introduce in Lcmma l5.1l to claiming that F{iy, R+) is independent of v. As 
we show in ProDOsition l5.6l this is indeed equivalent to the desired result. While we 
cannot rule out the possibility that this fact may be proved in some fashion trivial 
enough to be neglected, we have been unable to find a direct proof.^ 

Where we have made good use of a generalized version of the constructions from 
|CMM95] , is in deriving our Theorem 13.41 a result that has no parallel in their 
setting, but can only arise in the case of nonzero internal killing. If the killing rate 
is faster at infinity than near the origin, then escape to infinity is excluded: any 
mass that moves away from the origin will become relatively minuscule over time. 
This allows us to determine, when the limit of the killing rate at infinity is greater 
than the quasistationary killing rate, that the process does not run off to infinity. 

2. Definitions and assumptions 

2.1. The process. We consider a one-dimensional diffusion on the interval (0,r), 
where r is finite or infinite. We will always assume that is a regular boundary. 
If the right boundary is oo, we assume it to be inaccessible (that is, natural or 



As we discuss in the context of the proofs, some details of P. Mandl's paper needed to be 
fleshed out, to show convincingly that the methods could be generalized. In addition, some 
fundamental results about the generators of diffusions with killing had not been written down. 

"''The authors of ICMM95| have provided an alternative proof of their result, in a private 
communication, for their special case where killing occurs only at 0. Their proof of the lemma 
amounts essentially to a direct proof of their theorem, by proving Mandl's extra condition I34i . 
Indirectly, this also proves their Lemma 2. There is, as yet, no direct proof of their Lemma 2, 
hence no reason to think that the method proposed in ICMM951 can be made to work, even in 
the case of no internal killing. 
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entrance); if it is finite, we assume it to be regular. When r is not explicitly 
mentioned, it is assumed to be oo. 

The diffusion rate a is assumed continuously twice differentiable. The drift 
h G C^((0,r)) is assumed continuously differentiable. The killing rate, denoted k, 
is only required to be continuous and nonnegative. We write t^ for the first time 
the process hits state x. The time at which the process is killed, transitioning to 
the cemetery state d, is thus represented by tq. 

We may simplify the problem somewhat by assuming that a is identically 1. 
There is no loss of generality since (following |Fel52p we can replace Xt by Yt = 
F{Xt), where F{x) := / du/a(u), where xq is an arbitrary point in (0,r). The 
killing rate for Yt becomes K{F~^{Yt)), while the drift may be computed by Ito's 
formula to be 

(1) a(f-i(yo)" ^ ^'^^- 

Since a is positive on the interior of the interval, this transformation is finite. From 
now on, we will always assume, unless otherwise stated, that the diffusion rates of 
our one-dimensional diffusions are identically 1. We note that in this normalization 
the speed measure has density e^^^-' with respect to Lebesgue measure, where 
B{x) := 2 J^ b{z)dz, while the scale function is s{z) = J^ e~^'^^^dx. We denote by 
/3 the measure on M+ that has density e~^^^' with respect to Lebesgue measure. 

The diffusion semigroup acts on the Banach space C([0,r]) of functions contin- 
uous on the closed interval [0, r], with the supremum norm. (When r — oo this is 
the extended half-line. We will use the notations [0,cxd) and R+ interchangeably 
for the half- line including 0.) As a strongly continuous contraction semigroup, it 
has an infinitesimal generator, which is the closure of the operator 

(2) L4) = ]-(t>"{x) + h{x)(t>'{x) - n{x)4){x), 

acting on a subspace of the twice differentiable functions whose image under the 
operator is also in C([0,r]). (The subspace is determined by lateral conditions at 
the endpoints. Each regular endpoint requires a boundary condition to restrict the 
generator; the boundary conditions restrict the domain of the pregenerator, which 
is then closed to make a generator.) The dual space is a subspace of 'BV([0,r]), 
the signed measures with bounded variation on the closed interval [0, r], with total- 
variation norm. There is a natural inclusion of L^ into 'BV((0,r)), taking / g 
L^ to the signed measure with density / on (0,r). This identifies L^ with the 
subspace of absolutely continuous signed measures. Since the range of the resolvent 
is contained in this subspace, the semigroup may be restricted to L^. (This is proved 
in Theorem 13.3 of |Fel52| . which we extend to diffusions with internal killing in 
Lemma l7.ll Because the space 'BV([0,r]) includes masses at the endpoints, in 
general these could be retained in the image. In fact, though, when the boundaries 
are natural, entrance, or regular without "stickiness" , point masses at the endpoints 
are excluded from the image of the resolvent.) The adjoint generator is the closure 
of the operator 

(3) L*cf> = i</>"(x) - {h{x)c^{x))' - n{x)c^{x), 
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acting on the space of twice differentiable functions whose image is in L^, with 
boundary conditions imphed by the nature of the endpoints. (This description of 
the generator and the adjoint is given by W. FeUer in jFel52| .) 

The generator describes the behavior of the process in the interior of the state 
space. The behavior at the endpoints is determined by a boundary condition on 
the domain. If is a regular boundary, we need a boundary condition to describe 
the behavior. We assume that there is no "holding" or "stickiness" at the boundary 
— that is, the speed measure m satisfies m({0}) = — and we assume that b is 
continuous at a regular boundary. For this case. Theorem 19.3 of |Fel52| gives the 
boundary condition on the domain of £* as 

(4) {l-po)^{x)^po{^c^'{0)-b{0)cf,{0)), 

where < po ^ 1 is an arbitrary parameter. (Note: Everywhere that we have 
i(/)' — 60 for the boundary conditions, the sources we cite have twice that. Of 
course, this is arbitrary, since the normalization of pq is arbitrary; but it seems 
appropriate, because the sources are normalized with a^ = 2, whereas we have 
cr^ = 1.) It is a well-known fact that (1 — po)/po represents the rate of absorption 
at 0. If r is a regular boundary as well, the analogous condition holds there: 

(5) (1 - pr) lim e-^^^'Uix) = ~Pr hm (l<l>'{x) - 6(a;)0(a;)) . 
If r is an entrance boundary, there is the simpler condition 

(6) \m-,{U\x)-h{x)4>{x))={). 

a;— ►r Z 

(These boundary conditions are given in Fel52 for the forward semigroup and the 
adjoint, in the absence of internal killing. In ^IM65, the same conditions are derived 
for general killing, but only for the forward semigroup. We transfer this result to 
the adjoint in Lemma l7.1l ) 

The adjoint semigroup is not discussed in |IM65j . The arguments of |Fel52| 
carry over directly. We recapitulate them here, because an exact description of the 
adjoint generator — and of its domain, in particular — will be indispensable to the 
proof of Theorem l3.1l We define C to be the set of twice continuously differentiable 
functions on M+ satisfying the boundary condition Q at 0, and the appropriate 
analogue condition at r in the case of two regular endpoints. Given an operator 
A and a strongly continuous contraction semigroup Tt, we will write "^ generates 
Tt" or "Tt is generated by A" as shorthand for "A is closable, and the closure of A 
is the generator of Tt" . 

There are several additional conditions that we will need to impose on our dif- 
fusions. To begin with, we need 

(LP) £,* is in the limit-point case at oo. 

(For an account of the theory of singular Sturm-Liouville problems, and the concepts 
of "limit-point" and "limit-circle" in particular, see either chapter 9 of CL55 or 
chapter 5 of Yos91 .) Our main results will all assume condition (LP). This is fairly 
technical, but there is a simple condition on the diffusion parameters which suffices. 

Lemma 2.1. The condition 

(LP') liminfz"2(fe(z)2 + 6'(z) + 2K(z)) > -oo or k ee 0. 
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implies (LP). 

Proof. When k = 0, the proof is given in jManBlj . The limit-point property follows 
directly from the fact that 0(z) = e-^'^^ J^ e~^(^^da; is a solution to iL*0 = 0. The 
assumption that cx) is inaccessible implies that this is not in i^(/3). 
We make the standard transformation 

£0(.T):=e-^(-)/2i:(e^(-)/V(2:)) 

= ^n^) - libixf + b'{x) + 2^{x))^{x) 

The eigenfunctions of L may be transformed into those of £* simply by multiplying 
by e^'^'/^. Thus t, which is self-adjoint with respect to Lebesgue measure, is in 
the limit-point case at oo if and only if this is true of £*, which is self-adjoint with 
respect to e^'^^^dx. But Corollary 9.2.2 of |CL55j . together with the assumption 
(LP'), implies that L* is in the limit-point case at c». D 

Note that (LP') is not an onerous condition: it can fail only if —b'{z) becomes 
arbitrarily large (as z — > oo), while 5(z)^ remains relatively small. 

2.2. Eigenfunctions of the generators. For any A, we define (l)\ to be the unique 
solution to the initial value problem 

(7) i:*0A = -A0A with (/.a(O) = po and ^^'^(0) - 6(O)0a(O) = I - Pq. 
We define 

(8) A = supjA : 4>\{x) is nonnegative for all a;}. 

Similarly, we define ipx to be the unique nonzero solution to the initial value problem 

(9) L^x = -Mx with i^a(0) - po and ^^1(0) - 1 -Po- 

Note that i^xix) = e^^'^^^(j)x{x). If -0 is any A-invariant function — that is, a 
function ip such that Pti/j — e^^'^tp — then -0 is a multiple of 0a- The converse 
need not be true: that is, 0a need not be A-invariant. We follow Mandl in defining 
Xg to be the smallest point of increase of the spectral measure g for the operator 
£*. (We also follow Mandl in reversing the sign on the spectral measure; where 
the spectrum should be concentrated on (— cxd, — A^] — which would correspond to 
considering eigenvalues 0a with eigenvalue A, where A is negative — we have let 0a 
have eigenvalue —A, and put g on [Ag,(X)).) 

We follow Mandl in defining ^ to be the unitary map from L'^{(3) to L'^{g) defined 
by the eigenfunctions. That is, for / € L'^{(3), 

(10) */(A) = lim / f{x)<l>x{x)e-^^'''^dx, 

where the limit is understood in the L'^{0) sense. The existence of this limit, and 
the fact that it is unitary, are shown in Theorem 3.1 of CL55 . The measure (3 does 
not appear there, but the essential assumption is that the differential operator is 
self-adjoint. Since our differential operator L* (acting on the appropriate domain 
I), defined in sectionQJ is self-adjoint with respect to /?, the result carries over. 

An essential step in Mandl's method is his Lemma 2, the fact that A = \g. This 
depends only on the assumption that the operator is in the limit-point case at oo. 
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The proof is sketchy, though, in the original, which might leave one in some doubt 
of whether it can properly be extended to our more general setting. 

Lemma 2.2. Suppose (LP) holds. Then X = Xg is finite, and the set of X such 
that <j)\ does not change sign is (—00, A]. Similarly, the set of X such that ijj\ does 
not change sign is (—00, A]. 

Proof. Let I be the set of A such that the boundary value problem lO on [0, 00) 
has a nonnegative (and nontrivial) solution. Theorem 8.2.1 of |CL55| tells us that 
the r-truncated boundary-value problem has solutions for infinitely many A. By 
the r-truncated problem we mean the initial-value problem on [0, r] with the same 
condition at 0, and the additional condition (f){r) = 0. These A may be ordered 
as > — Ao > — Ai > • • • , where lim„^oo A„ = 00, and the eigenfunction c/jx^ has 
exactly n zeroes on (0,6). Any A for which there is a solution to the r-truncated 
problem for some r is automatically not in 2. This implies that the complement of 
X is unbounded above. 

Suppose now that we had three values A+ < A* < A_, with A+,A_ G X, but 
A* ^ X. Since the zeroes in any solution to the differential equation are discrete, we 
may define an extended real-valued function C(A) on [A+,A_] as the smallest zero 
of (f>\, or 00 if the function has no zero. On any compact set of z and A, (j)x{z) is 
uniformly continuous in z and A (by Theorem 1.7.5 of iCL55j V This implies that C 
is continuous. Since C('^+) = C(-^-) = +00, and since C(-^*) is finite, there must be 
A' < A", both in [A+,A_], such that C(A') = C(-^")- But this would imply that the 
6-truncated problem for b — C(A') has a nonnegative solution for two values of A, 
namely A' and A". This contradicts the above-mentioned Theorem 8.2.1 of |CL55| . 
The same argument shows that we cannot have A+, A_ ^ X, but A* € X. 

We may conclude that X is either empty or an interval unbounded below. Mandl 
shows clearly that Xg € X. The proof depends only on the fact that the initial- 
value problem is in the limit-point case at 00, which we have assumed. On the 
other hand, Xg is a point of increase for g defined as a limit as r ^ 00 of step 
functions with jumps at the eigenvalues for the r-truncated problem. There must 
be A > Ag arbitrarily close, such that A is an eigenvalue for the r-truncated problem 
for some r. These A are not in X. Thus, X = (— cxd, Xg] = (—00, A]. 

The same argument works for the boundary- value problem ^. D 

2.3. Asymptotic properties. We say that the process Xt has an asymptotic 
killing rate rj when started from the probability distribution v if 

(11) lim P^jra >t-|-s|ra >t| =e"''^ 

for any positive s. A trivial consequence is that 

(12) 77 = - lim f"4ogP^{ra > t|. 

If G Co([0, cxd)) is nonnegative, we say that the process Xt converges from the 
initial distribution v to the quasistationary distribution </> on compacta if for any 
positive z, and any Borel A C [0, z], 

(13) Yiu,^,{Xt^A\Xt<z]^ \Y^ . 

t^oo ^ j^ (j){z)dz 

It should be noted that this property, while interesting, is far short of what one 
would like to prove. In particular, it leaves open the possibility that most of the 
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mass — an amount which falls off at a rate strictly slower than e~-* — runs off 
to oo, even while the part of the distribution which stays below z, for any fixed z, 
converges to the distribution with density (j)\ . We say that Xt converges from the 
initial distribution v to the quasistationary distribution (p if H13|l holds with z = oo; 
that is, for any Borel subset A C [0,oo), 

I , [a (t>(z)dz 

(14) lim P,{Xf G A ra > f} = -^^rjT^TT- 
t->oo Jo V\z)d-z 

We say that Xt escapes to infinity from u if 

(15) lim P^jXt < zlra > t| =0 

for all z G M+. 

An elementary consequence of the definitions is: 

Proposition 2.3. If Xt started at v converges to the quasistationary distribution 
4> on compacta, and J^ 4>{y)dy = oo then Xt started at v escapes to oo. 

Proof. Given any z' > z > 0, 

limsupPj.{Xt <z\To>t}< limsupPj.{X( <z\Xt<z'} 

t—^oo i— >oo 

^ lo Hy)dy 
lo 4'iy)dy 

Sending z' to oo proves the result. D 

2.4. The initial distribution. The convergence properties may depend on the 
initial distribution. In particular, if the tail of the initial distribution is too heavy, 
there may be mass wandering gradually in from the tails more slowly than the 
decay rate of the mass starting on any compact set. (In the case of countable-state 
discrete Markov chains, conditions on the initial state leading to convergence to 
the quasistationary distribution are discussed in |SVJ66j .l We adopt the following 
sufficient condition from Mandl: 

(ID) The initial distribution has a density / G L'^iP) with liminf ^/(A) > — oo. 

(Unless otherwise indicated, densities are always meant to be taken with respect 
to Lebesgue measure.) Since we are concerned with asymptotic properties, there is 
no change if we replace Xt by Xt+s, for some positive s. Thus, we may substitute 

If Xq has distribution i^, then 3s > for which the distribution of Xg 
(^°') has density / G L^(/3), with liminf ^/(A) > -oo. 

A], A^ 

We will use (IDC) to denote the condition that a distribution satisfies (ID') and 
has compact support. The following is a trivial consequence of the definition Hl()|l 
of ^: 

Proposition 2.4. // the initial distribution has a density with compact support, 
then condition (ID) is satisfied. 

It seems plausible that condition (ID') should be satisfied for any initial distribu- 
tion with compact support (which seems to be implied by comments in |Man61j ). 
but we have not been able to show this. 
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2.5. The Lo function. We define the spaces 071^ (M+), comprising probability distri- 
butions on [0,00) with compact support, and 37l<i(R+), comprising subprobability 
distributions on E+. When po = (that is, in the case of complete killing at 0), 
we exclude the point Sq from 9J15(R+). We supply 9Jl<i(R+) with the topology 
of weak convergence on compact sets, while 3Jt^(R+) has the stronger topology in 
which Un ^> u ii and only if the sequence converges weakly, and IJ^^ Supp Vn is 
bounded. (That is, the open sets are generated by weak open sets, together with 
{v : Supp V C (0,i?)} f or R > 1.) 

We imitate |CMM95] in defining functions 

(16) ^t(^):- p^|^^^'| for.e9Jl-(R+), 

mapping compactly supported probability distributions on M+ to R+. We also 
define 

(17) Ftii^,-):=P,{Xte-\Ta>t}. 
When 1/ — 5x, we will also write Ft{x, A) and (jJt(x). Note that 



(18) tJt(i/) = / LUt{x)diy{x) 

Jo 

for any compactly supported i^. Since uJtii^) is nonzero (except in the case of 
complete killing at 0, when v = Sq), we also have, for all i' 7^ So, 

, ■, , s L Ft{x,-)u)tix)dv(x) 

Jo ^t{x)dv{x) 

We note that a version of ujt (where the argument was confined to point masses) 
was introduced in ;CMM95) . 
Finally, we define 

(20) at(i/,r) i^P^ra >t + r|ra>i} = Ft{v,dv)fy{TQ > r] . 

For some results, we will need the following condition: 

(GB) Vs > 0, E:, s\rpujt{Xs) < 00 

t>o 

This seems to be an abstract condition, which might be hard to verify. In fact, 
though, it is a simple consequence of some fairly weak conditions on the diffusion 
parameters b and k. 

Consider first, the situation that k is monotonically increasing and there is com- 
plete reflection at (as in our example of sectional. We can couple a copy of the 
process Xt started at 1, to a copy Xj. started at a; > 1 in such a way that Xt < X[ 
while both processes are alive. At any time when both processes are alive, the 
killing rate must be larger for X' than for X, so we can also arrange the coupling 
so that X' is killed before X. This implies that iOt{x) < 1 for all x and all t, and 
hence condition (GB) is trivially satisfied. 

If K is not monotone, we need to restrict its growth, as well as the growth of the 
negative drift. We note that there is a general bound on uJt{x)/uJt{y), obtained by 
factoring Py{Td > t} over the time (if any) when the process started at y first hits 
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<:. We have 



This yields 
(21) 

(22) 



iOt[x) _ P^lrg > t} ^ PxJTa > t) 



^*{y,x) < — -r^ <uj^,(x,y), and 



UJtix) 



f^tiy) 



^t{y) 

< w*,(x, y) := uj^{x, y) ~ uj^{y, x)~^. 



For ^ we use the fact that ii A < B < C mid A < 1 < C, then \B - 1\ < C - A. 

Lemma 2.5. Suppose there are positive constants n^, 6,, &,*, anrf (3 such that 
either 

(GB') \b{y)\ — b*y ^''^'^ '*(y) — '^*y fory sufficiently large 

or 

(GB") 
— 6*y < h{y) < b^y, b'{y) < b^^,y and n{y) < K^y for y sufficiently large. 

Then for all s,t,x £ M+, 

(23) E,supcJt(X,) < E,w,(X„l) < c». 

t>o 

Proof. Integrating by parts, we need to show that 



Hm 



^x{Xs > W} 

'^ Pi{r„ < Ta} 



"JX, > z} fd ^ , , , , 



\{tz < TdY \dz 



.{>■) 



is finite. For A nonnegative, define Tq to be the minimum of rg and an independent 
exponential stopping time with expectation 1/A. Then 

Px{^. >z]< P^JT, < 4 < e^^P^r, < T^^)}. 

Since ri ' < rg, it suffices to show 



(24) 



"x{t. < r^^)} (d ^ . ^ , ^ 



i{tz < Ta}2 \dz 

Observe now that Py {tz < Tq } for y e (0, z) is the Laplace transform of Tz . On 
this interval it coincides with eigenfunction ijj-x{y)/'4'-xiz) oiL. What we need to 
show, then, is that 



(25) 

for some positive value of A. 
We define 



(26) 

Note that g\ satisfies 
(27) 



dz < oo 



9\{x) 



tp-xiz) 



i^-x{x)' 



gx{x) = ^m (Aa;)-ip,_A4^- > ^a }' 
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from which it foUows that g\{x) is positive, and nondecreasing in A. With this 
notation, the condition (|25|l becomes 



(28) j go(^)expM (.go(y)-gA(2/))|d2< 



CX). 



We now need an upper bound on g\{x). Fix any positive 5. We start the diffusion 
at a; — e, where e is assumed smaher than 5. Note that if we increase the kiUing 
rate on (x — 5,x) to k := supy^r^_g,^\ K{y), and kill the process when it reaches 
X — 6, that increases g\{x), since the probability of being killed before reaching 
X is increased. If we then decrease the drift rate on (a; — d, x) to the constant 
b :— irifye[x-s,x) f^iy)! that further increases g\{x). We can then apply equation 
2.3.0.5 of BS96_ to compute 

g^(x) < ^ ( ^^b-^JlM^^f^^] < -b+acoth{ad) < -2(bA0)+V2k + 2X+6-^ 
ax \ sinh(cH5) J x=o 



_2 2 

where a ~ \j b +2k + 2X. If (GB') holds, we fix, arbitrarily, (5=1, and we see 
that for x sufficiently large 

(29) gx{x) < 2b^x + y/2^x + V2\ + I. 

If (GB") holds, we take 6 = 1/x. Then for x sufficiently large 

(30) gxix) < -2{b{x) A 0) + 26,,x + ^/2^a; + V2X + x. 

Define h{y) := gx{y)~go{y), and k{y) := gx{y)+go{y) + 2b{y). Using either ^ 
or H30|l . we see that 



(31) k{y)<Cy + V2X + C', 

where C and C" are constants independent of A. 
From the equation ^ we derive the relation 

(32) h'{y) = 2\~k{y)h{y). 
This equation has the general solution 

(33) h{z) = h{l)e- ^1 '=(^)'^y + 2A / e" ^y '^^^^'^^dy. 

Let C" = max{V2A + C", 1}. By ^, for z > 1, 

h{z) >2X [ exp|-C(z2 - j/2)/2 _ C"{z ~ y)]dy 

Jz~l/C"z *■ -* 

>iV2X-C')e-^^'-'z-\ 
Thus, 

exp {^' (goiy) ~ gx{y)) } < z-(^-^>-^^-^ . 

Since go grows at a polynomial rate, this proves H28() for A large enough. D 
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3. Outline of the paper and statement of main results 

We will consider here mainly the case where is a regular endpoint, and r — oo 
is an inaccessible endpoint. The case of two regular endpoints is an application of 
basic results about compact operators, and differs only slightly from the standard 
K = case, such as is presented in |Pin85j . We present the more general proof 
only briefly. The case of two inaccessible endpoints also would offer no additional 
novelty, except to complicate the notation. 

Mandl's argument for diffusions on R+ proceeds in two steps. In the first step, 
he proves that Xt converges to the quasistationary distribution 0a on compacta. 
In order to remove the restriction to compact sets, Mandl imposes a condition 
equivalent to 

(34) / e^(^)dz < oo and lim / I / (f>x{y)dy\ g{dX) = 0, 

implying a strong drift away from oo. Collet et al. |CMM95] offer a different 
approach to removing the restriction to finite z. Under the assumption that the 
process is killed instantaneously at 0, and still with k identically 0, they claim to 
show that they can dispense with the assumption H34I) . Further results for this case 
are developed in [MMOlj . 

We wish to extend these results to nearly general killing, and more general bound- 
ary conditions. For the first part, the difhculties are merely technical. Mandl's 
proofs rely on W. Feller's careful analytical treatment of the diffusion semigroup 
|Fel52| . The relevant parts of Feller's work assume k = 0, and the corresponding 
results of Ito and McKean |IM65| . which do allow for general killing, are not suffi- 
ciently elaborated to be applied in this context. One part of our work will thus be 
to extend Feller's results rigorously to the case of general killing, and then to show 
that Mandl's methods still work in this setting. This also allows us to dispense 
with Mandl's requirement that the inaccessible boundary be natural; at the same 
time, we are forced to impose the minor additional condition (LP). Our first result 
is then 

Theorem 3.1. If condition (LP) holds, and the initial distribution satisfies (ID'), 
then the process Xt converges to the quasistationary distribution (j)\ on compacta. 

These technical extensions also allow the case of two regular boundaries with 
internal killing to be handled. We define A'-'^-' as the 

A^*^' = infJA : ipx satisfies condition (0)}. 

Theorem 3.2. Suppose that and r < oo are both regular boundaries. If both 
boundaries are reflecting and k = almost everywhere then A = 0; otherwise, A is 
positive. For any Borel subset A C (0,r), and any starting distribution jjl, 

P5) ■i..A,,„{.Y.C.4} ^''°'>'tt;':'-{f"'''''^ . 

t^oo jp e ^(y)(t){{y)dy 

We note that the assumption that b is continuous at the boundary could be 
dispensed with in the case of two regular endpoints. In this case, though, we can 
no longer appeal to general principles (the Krein-Rutman Theorem) to show that 
there is a positive eigenfunction at the bottom of the spectrum. We would then 
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need to impose this as a condition: 

(36) 4>\{r) is nonnegative. 

When the interval is infinite, the second part of the problem, proving convergence 
without conditioning on a compact set, requires a new approach. We split this part 
again in two. One part, which generalizes the results of |CMM95] . is 

Theorem 3.3. Suppose (LP) holds, and either 

liminf k(2) > \ or limsupK(z) < A. 

Then either Xt converges to the quasistationary distribution 4>x from every com- 
pactly supported initial distribution, or Xt escapes to infinity from every compactly 
supported initial distribution. In the case liminf 2,^00 k(-z) > A; ^t converges to the 
quasistationary distribution (px if and only if J^ 4'\{y)dy is finite. 

Note that we have already shown, in Proposition 12.31 that the process must 
escape to infinity if J„ <p\{y)dy — cxi. When the integral is finite, though, there 
may be ambiguity. 

We remark here, as well, that Lemma l^l^ tells us that the process has asymptotic 
killing rate A when it converges to (/>>• li K := lim^^oo k{z) exists and is finite, 
the process has asymptotic killing rate K when Xt escapes to infinity. (Unless 
otherwise indicated, when we say that liniz^oo k{z) exists, we always mean to 
include the possibility that limz^oo k{z) — 00.) 

The general idea here is that the process cannot split, with some mass going off 
to infinity and other mass staying near the origin, if the killing rates differ. If the 
killing rate converges to a limit which is larger than A, we can be sure which side of 
the dichotomy we are on, namely, convergence to the quasistationary distribution. 

Theorem 3.4. Suppose the endpoint at cxd is natural, conditions (LP) and (GB) 
hold, and K := liniz^oo k(z) exists. 
(i) // K 7^ A, then either 

• Xt has an asymptotic killing rate rj ^ X and converges to the quasi- 
stationary distribution 4>\ for every compactly supported initial distri- 
bution; or 

• Xt has an asymptotic killing rate rj — K and escapes to infinity for 
every compactly supported initial distribution. 

(ii) If K > X, then rj ^ X. 

(iii) If K = X, then Xt has an asymptotic killing rate rj ~ X^ K . 
(iv) Regardless of the relation between K and X, we have 

(37) lim .,(.) = hm M^^ = #1. 

Note that this result dispenses with the constraint (ID') on the initial condition. 

In section El we present, as an example, a continuous analogue of the Le Bras 
process. The Le Bras process, which was mentioned in section^ is simply a Markov 
process on the positive integers, in which the transitions are only single steps up, 
with the transition rate proportional to the state. The process is killed at a rate 
which is also proportional to the current state. 

For the sake of clarity of exposition, we have placed all lemmas whose proofs are 
merely technical, and do not shed light on the central issues, in section [3 at the 
end. 



Generator 


Space Norm 


n* 


L' ll/lli = /o"'l/WM-- 


Q 


L'iP) ii/ii2 = (/o"i/(z)mdz))'^' 


n 


i = L^nL^P) ||/||dl/l|i + ll/l|2 
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4. Convergence on compacta 

4.1. Proof of Theorem 13. IL This result generalizes Theorem 2 of jManBlj . The 
proof there assumed k = 0. We follow the general lines of Mandl's proof, substi- 
tuting stronger (or clearer) arguments as needed. 

We begin by defining three differential operators, il* , fl, and il. All three act in 
the same way, as L* , on different subsets of C, the set of continuously twice differ- 
entiable functions on R+ satisfying the boundary condition Q) at 0. The subsets 
are assigned different Banach-space structures. We use the following notation: 

Domain 

■D = {/e enLi : nf eL^} 

7) = {fjE en L^iP) -^nf E L^m} 

■D = {fee:Qfee} 

We note that all three domains are dense in C: any function in C may be approx- 
imated in any of these three norms by a function that vanishes, along with all its 
derivatives, outside of a compact interval. Mandl's Theorem 1, which we need to 
prove, states that all three of these operators generate strongly continuous contrac- 
tion semigroups, and that the restrictions of all three to i are identical. 

The operator fl* is precisely the one analyzed in Lemma l7. II There we showed 
that it generates the adjoint semigroup of our killed diffusion on L^. 

The content of Problem 9.12 of |CL55| is that for any a not in the spectrum of 
ft (in particular, for any positive a), for / G i^(/3), 

m ,„-.,-.;. *-.(iM)^ 

(Note again that the sign on A is reversed, because we have flipped g onto the 
positive half-line; in particular, we have defined (j)\ to have eigenvalue —A.) Then, 
since ^ is unitary. 



a\Ua - n)-' f\\ = 



a + A 



*/ 



<ll*/ll = 11/11 



(The inequality depends on the fact that the spectrum of f2 includes no positive 
values.) This means that Cl is dissipative, proving (by the Hille-Yosida Theorem) the 
existence of a strongly continuous contraction semigroup on L^ (/3) , whose generator 
is an extension of ft. To prove that the generator is actually the closure of Cl, we 
need to show (by Theorem 2.12 of |EK86| ) that the image of T) under a — 51 is dense 
in I). Let S be the subset of L'^{g) consisting of functions with compact support. 
Consider {a — (l)^^'^^^{S). All functions in this image are bounded integrals of 
4'ix, A) (with respect to A), hence they are twice differentiable, and the continuity 
of the diffusion coefficients implies that they are continuously twice differentiable. 
Since ^^^(5) is dense in L^(/3), this proves that D is a core. 

Now choose f G I, and any A > 0. Then (by Proposition 1.2.1 of [EKSSp there 
exist unique solutions ui € L^ and U2 & L'^iP) to 

(39) Xui - n*ui = /, 
and 

(40) Am2 - nu2 = /, 
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with A||wi||i < ll/lli and A||u2|l2 < H/lh- (Here || • ||i is the L^ Lebesgue norm, and 
II • II2 is the L^ norm with respect to /?.) We show now that ui = U2- 

It suffices to assume that / is nonnegative, since we may decompose / into 
positive and negative parts. Since ui = /^ G*^{x,y)f{y)dt, and G*^{x,y) > 0, the 
function ui is nonnegative. Equation H38|) implies that U2 is nonnegative as welL 

As we have aheady shown in the proof of Lemma l7.4l there exists a nonnegative 
solution u to the equation (L — k)u — Au = 0, satisfying the boundary condition 

The function v{x) := e^^^'g{x) is then a nonnegative solution to {H* —k)v — \v — 0, 
satisfying the boundary condition 10} . Since the solutions to ifS^ and ff^ are 
unique, it follows that v is neither in L^ nor in LP'{j3). By uniqueness of solutions 
to ordinary differential equations, we must have ui — U2 — cv, for some c. If c > 0, 
then ui > cv > 0, so wi would not be in L^. On the other hand, if c < 0, then U2 
would not be in L'^[(3). Consequently, c = 0, and ui = U2 is in L Furthermore, 

||u,||, = ||ti,||l + hd|2<||/||i+ 11/112=11/11,, 

SO f2 is dissipative. It follows from the Hille-Yosida Theorem |EK86I Theorem 1.2.6] 
that Q. generates a strongly continuous contraction semigroup, for which £ is a core. 
At the same time, we have shown that the resolvents of all three generators are the 
same on £, so the restrictions of the semigroups to t must also be identical. 

We note, now, that H38|l implies, by identity of resolvents, that the semigroup 
Pt generated by tl maps under ^ onto the semigroup Qt on L-'{q) defined by 
QtgW = e~^*g{X). Consequently, for f € £, 

(41) *P;/ = ^Ptf = e-^**/. 

We may assume, without loss of generality, that (ID) is satisfied, rather than 
(ID'). Let / be the initial density. For any bounded Borel set A, by H41|l and the 
unitarity of ^, we have 

(42) F{Xt eA}^ J lA{x)P:f{x)dx = j e-^'UWg{X)dg{X), 
where g :— \I'/ and 

UiX) := ^(e^lA)(A) = / e^(^VA(2:)dx. 

JA 

Note that 'Ca(A) is an integral over a finite interval, hence is continuous. For any 
h G L^{b) define (following Mandl) 

I{t,h) := / e-^'h{X)dg{X). 

Following 142|l . it suffices to show, for bounded Borel sets A and A' , 

(43) lim ^(^'g^-^)^^(^). 

For h : (A, oo) -^ M+ and A* > A, define I{t,h,X^) :— X(i, /i1[a,a.])- We note 
first that for any ft-i, ft.2 G L^{q), with /i2 positive, and oo > Ai, A'^, A2 > A, 



I{t,hi,Xi) I{t,hi,X[) 



I{t,h2,X2) I{t,h2,X2) 



< e 



(A-AiAA, AA2)t -^A 



4°°|fei(A)|dg(A) 
j;^ h2iX)dgiX) ' 
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SO that 

lim sup :7- — ; — ^ and lim mi — - — — — ■ 

are independent of Ai and A2. Moreover, we have the bounds 

(44) hm sup , ,,, > hm sup — - — ; -— > hm mf — - — ; -— > hm inf -p— -— 

AiA /12(A) - t-,00 X(i,/l2,A2) - t-00 X(i,/i2,A2) - AiA /12(A) 

In particular, assumption (ID) imphes 

(45) sup lim sup • < 00, 

A,>A t^oo -i^r, i,A*J 

where g" = max{— g, 0} and g^ — max{g, 0}. 

Suppose that / has compact support. Then g is continuous, and 

/•oo 

5(A) - / I{x)4>x{x)dx > 0. 

"'0 

Thus, for any finite A* > A, 

lim|^^^=eA(A)ff(A)>0, 
t^oo _L(t, 1, A,) 

from which l|43|l is immediate. 

Now, consider a general density /. Choose zq such that Jq° f{y)dy > 0, let 
/o = /l[o,2o]' ^^^ let 50 = */o- Then 

(46) hm inf ^^ > lim ^^^iMa). > Q. 
^ ^ i^oc X(i, 1,A*) ~ *^°°2:(/, 1,A,) 

We claim that 

(47) 5 := 1 - sup lim sup jfj^'^ } > 0. 

Suppose this is true. Then 

,. At,\g\) ^,. f I{t,g+,X,)-I{t,g-,X ,)Y^ 

hmsup — r— < hmsup _ 

t^oo I{t,g,K) t~,oo V 2:(t,5+)+2(t,5 ) 

(48) ^ f. .. I{t,g-,X. r -' 

< 1- hmsup 

V t^oo I{t,g+) 

<s-\ 

We may then conclude that for any continuous function h : [A, oo) -^ M, 

h{X) - 5-^ sup |/i(A) - /i(A)| < liminf J^^'^y, 
Ae[A,A.] *^°° I[t,g,X^) 

I{t,gh) 

- limsup— — - 

t^oo l[t,g,X^) 

<h{X) + S-^ sup \h{X)-h{X)\ 

Ae[A,A.] 
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Thus, for any A* > A, and continuous hi and /i2 with /i2(A) > 6^^ sup;^gr^ ^ ^ |/i2(A)— 
h2{X)\, 

ShljX) - SUP;,g[j, ;,_^] |fei(A) - fel(A)| 
J/l2(A) +SUP;,g[2,,A.] I^2(A) - /l2(A)| 

^ y ■ r^it^ghi) i{t-,ghi) 

< um mi — < Imi sup — - 

t^oa I{t,gh2) t^oo I{t,gh2) 

^ (^fellA) + SUP;,g[^_^^] |/ii(A) - fei(A)| 

- Sh2{X) - sup;^g[^_;,^] |/i2(A) - /i2(A)| ■ 
Sending A* to A, we see that 

^ ' t-*ooJ(t,g/i2) /l2(A) 

Applying this with hi = ^a and /i2 — ^a' completes the proof of (021 • 

It remains to prove the claim (|T7| . By continuity of ^a, and l|^ . we may find 
A* > A and a positive constant e such that for all A* e (A, A*), 

liminf^^f4^1^-2 sup |a(A)-a(A)|- sup .g-(A) > elimsup^^*'''^ ^ 



t^oo X(i, 1,A*) Ae(A.A*) Ae[Aa.]' t^oo 2:(t, 1) 

Define (5(e) := e^A(A)/(l + e^AiX))- Suppose, for some A* < A*, there were a 
sequence t„ — > oo such that 

n^oo i,(t„,5+. A*) 
Then 

(50) hm sup < j^ hm sup . . and 

/riN ,. 2:(t„,.g) 5(e) I{tn,g~) 

(51) limsup— — ^ < r-r^limsup: 



n— ♦oo -i-(En,ljA*) 1 — d(e) n^oo -^(tnjljA*) 

We may conclude that 
hminf^/^-'^y <a(A)liminf: ^^'"'^^ 



n-.oo X(t„,l,A*) n^oo I(i„,l,A*) 



+ sup |^^(A)-^^(A)|.-^— hmsup;^^*"'-"' 



AG(A,A*) -L — ^lej n-+oo -i(,E„,i,A*j 

SO that 

. „ 2:(t„,5) . e ,. I{tn,g-) 

hm mf — -^ > hm sup ; 



n^oo I(i„,l,A,) ~ ^a(A) n-too X(i„,l,A*)' 

contradicting H51|l . We have thus proven the claim 14711 for A, sufficiently close to 
A. By monotonicity, this completes the proof. 

4.2. Proof of Theorem 13.21 Let Ga{x,y) be the a-Green's function for the op- 
erator £, with respect to Lebesgue measure. The Green's function for £* is then 
given by G%{x,y) = Ga{y,x). Equation 6 of section 4.11 of |IM65j gives us an ex- 
plicit representation of the Green's function Ga{x, y), from which we can derive (by 
Theorem VI.23 of |K,S72j ) that the operator on L^(/3) generated by the kernel G* 
is compact. It is also self-adjoint with respect to the measure /3. By the Hilbert- 
Schmidt Theorem (Theorem VI. 16 of |R,S72j l the spectrum is discrete, and has 
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only as accumulation point. Furthermore, there is a complete orthonormal basis 
of eigenfunctions. These are also eigenfunctions for the infinitesimal generator; the 
eigenvalue —A for the generator corresponds to —^ for G^. Let — Aq > — Ai > • • • 
be the eigenvalues, and 4>k the eigenfunction of the generator Q corresponding to 
Afc. Since the eigenfunctions are orthogonal with respect to the positive function 
g-B{x)^ no more than one may be nonnegative. Either by assumption, or (if the 
drift is finite at the boundaries) by the Krein-Rutman Theorem in the form derived 
as Theorem 2.1 of |Nus84| . there is a nonnegative eigenfunction, and it corresponds 
to the top eigenvalue — Ao = —A. (With a slight further restriction of the conditions, 
we could instead apply the more elementary result. Theorem 8.2.1 of jCL55| .) 
If h{x) is a smooth function on [0, r] with 



/i||2 := / h{xf(3(dx) <oo. 



we can write 






A:=0 ^"^ ^ fc=0 

(The sum need not converge pointwise, but in L'^{I3).) Let 

K 

/if(a;)=^e-^'=*afe0fc(x). 

fc=0 

If the initial distribution has a smooth density /i, then by the Cauchy- Schwartz 
inequality (following an equivalent computation on page 561 of |Man61| ). 



lim 

t — ^oo 



e^*P{Xt e A} - ao / (l>x{x)dx 



A 



hm lim / e-^(-)/2|e^*/if -ao0o|e^("'/'da; 

t— ►CXD K—*oo J j^ ' ' 

< lim lim II e^* /if -ao0o II • / e^^^^dx 



'0 

= 0. 

If the initial distribution does not have a smooth density, we start instead with 
the distribution at a small finite time e. 

5. Convergence on R+ 

5.1. The limit functions tu and F. 

Lemma 5.1. Given any increasing sequence i* — *■ cxd, we may find a subsequence 
t„ — > oo with the following properties: 

(i) There is a limit uj^^"'{x) = linin~>oo 0Jt„{x) in C(R^), with convergence 
in the topology of uniform convergence on compacta. The function to is 
nonzero everywhere, except when the process is killed instantaneously at 0, 
in which case w(0) = 0. 

(ii) There is a limit F^*") = lim„_oo i^t„ in C{^l{M+),m<i{^+)), supplied 
with the topology of uniform convergence on compact sets. 

When there is no need to emphasize the particular subsequence, we will drop the 
superscript (tn) from the notation. 
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Proof. We consider first tlic functions Ft and cut restricted to (5-measures. Tlie 
topology on 9Jlf (IR+) restricts tfien to the usual topology on R+. 
For integers n > 1, t e [0, oo), and a; e [n, n + 1], define 

We show that for each n the set of functions „Wt is uniformly equicontinuous. The 
bound (|21|l immediately tells us that 

nUJt{x) < P„{r„+i < ra} < oo 
for all X, so the functions are uniformly bounded, li x,y G [n, n + 1], then 

\nUJt{x) - „UJt{y)\ <Pn{Tn+l <To} Vx{Ty<To} VPy{T^<ra} -1 

Since the diffusion rate is constant, and the drift and killing are bounded on [n, n+l], 
this bound goes to zero uniformly as jx — y| — > 0. 

By Ascoli's Theorem |Edw95l Theorem 0.4.13], the set of functions {ni^t}t£R+ 
is relatively compact. Thus, any sequence t'^ —> oo has a subsequence (U), such 
that ni-^ti converges uniformly to a function w*^"-' : [71,71 + 1] -^ M+. Furthermore, 
we may take the sequence for nWt to be a refinement of that for n~i^t- By taking 
a diagonal subsequence, we then get a single sequence U such that n'jJti converges 
uniformly to „a; for every 71. By an identical argument, we can also assume that 
QUJt := fx{Td > t}/Pi{Ti) > t} for X G [0, 1] also converges along the sequence 
ti. (We treat the case a; < 1 separately only because killing at could make it 
inconvenient to place PoJTa > t} in the denominator.) If we now define cut by 
uJt{x) = Fx{tq > t}/Fi{Tg > t} for all x G [0,cx)), we see that ujt is a product of 
finitely many functions n^t, so that tu = linii^oo '^t; exists as a positive continuous 
function, though the convergence need only be uniform on compacta. 

We note, as well, that for all y e R+, by (EH), 

w(y)>F,y{ri<ra}, 

which is positive, except when y = and there is complete killing at 0. 

We continue with Ft. We note that it will suffice to prove (ii) for the measures 
nFt, defined by 

nFtix,A) ^Ftix,An[-,n]); 
n 

the measure F{x, ■) is then defined as the limit along a common refinement of the 

subsequences defining „F(a;,-). We need to show that the set of functions Ft is 

relatively compact. Since relative compactness is a topological property, it will 

suffice to prove equicontinuity in any convenient metric that metrizes the compact 

weak-* topology. We choose the Wasserstein metric |Rac92| . given by 

d{n,fi') =sup{|^(g)-^'(5)|}, 

where the supremum is taken over functions g : M+ -^ [0, 1], supported on [—,71], 
and with Lipschitz constant < 1. 

Since the space is bounded, we need only to show that for any x G [0,oo) (or 
X € (0,00), in the case of complete killing at 0), 

(53) lim snpd{nFt{x,-),r.Ft{x',-)) =0. 

x'^x t 
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Let Rx'.x be the (subprobability) distribution of Tx, when the process is started 
from x' . Let g be any appropriate test function. 
By the strong Markov property, for any toi 



^x'{to > t} ^ 
Px{rs > t} 



nFt{x,dy)g{y)- / nFt{x' ,dy)g{y) 

C* Rx'Mds)\Exg{Xt_,) - ExgiXt)\ 



< 



r{Td > t} 



+ i?x',x([io Ai,oo]) 



Ex' [g{Xt) Irg > tkTx > to At] -EJgiXt)\T > t] 



< LO^^:{x,x') + / Rx' ^xids)Px{Xt-s > n+ I and Xt < n\Td > t^ 



Rx'.x{ds)^x [(1 A \Xt - Xt-s\)l{Xt-s < n + 1} I ra > i 

+ Rx' ,x{[to At,(xi\)iEx'[l /\ \Xt-x'\ \Td>tkTx>toM] 

+ Ex[l^\Xt-x\] +|a 

Note that this bound is independent of g. It remains, then, only to show that the 
supremum over t converges to as a:' ^ a;. By regularity, Rx'^x converges to (Sq- 
Since lj*, is continuous, and (jj*,(x,a;) = 0, the first term goes to as x' ^ x. The 
second term is bounded by 

(54) / R^{ds) sup fy{Xs < n I Ta > s] . 

Jo V>n+1 

The third term is bounded by 

nt 

(55) / R,{ds) sup Ey[{lA\Xs~y\) sup Pj^jra > s}"\ 

Jo !/<"+l y<n+l 

The Bounded Convergence Theorem implies then — since X is a Feller process — 
that lf51|l and (|^ both converge to as i?jy ^ Jq- 
We may conclude, for any positive tg, that 



lim sup 

x' — >X I 



nFt{x,dy)g{y)~ nFt{x' ,dy)g{y) 



< lim supi?2;'.2;([to Ai, cx)]) 

x'~*x t 

X (e^, [1 A \Xt - x'\ I Ta A T^ > t] + E^l A \Xt ~ x\] + \x - x'\ 



< lim 2Rx',x([to,oo]) + lim sup Rx'.x([t,oo])Ex[l A\Xt - x\] 

X'^X x'^Xf^f.^ 

+ lim supi?a;' a;([t, oo])E2;'[1 A\Xt- x'\\to > tkTx > to At]. 
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By continuity of the diffusion parameters, we may find an increasing function h(t) € 
[0, f] such that limt^o h{t) — 0, and for all x' < x + 1, 

Ex> [l A \Xt - x'W Tg > t k T^ > to A t] < h{t) and 

E.,[lA\Xt^x\\TQ>t] <h{t). 

The bound may then be turned into 

lim sup2i?:r',x([i,oo])/i(i) < 2/i(ti) + 2 lim R^,J[ti,oo\) <2h{ti) 

for any ii > 0, which goes to as ti ^ 0. 

It remains to extend the result to all v G SfJIf (R+). We define u; and F by the 
formulae fH^ and fH^. Since u}{x) > for a; > 0, w(z^) is nonzero, so that F(i/, •) 
is well defined. Since uj is bounded on compact sets, for any sequence tn such 
that lim„^oo wt„ (a;) = w(a;), we also have Vnnn^oo ^t„{^) = '^i'^): by the Bounded 
Convergence Theorem. If ;/„ is a sequence of measures converging to v in 971J(R+), 
they converge weakly and have their supports contained in a common compact set. 
Since uj is continuous, it follows again from the definition of weak convergence that 
lim„^oo / w(a;)(ii/„(a;) = J uj{x)diy{x). 

Now, consider any sequence i„ -^ oo, such that lim„^oo Ft^{x, •) = F{x, •). The 
same bounded convergence argument shows that lim„_»oo Ft^ (v, g) = F{v, g) for 
any v e 9Jt^(M+) and any continuous g : R+ — > [0, 1]. Suppose Vn is a sequence of 
measures converging to v in 97lf (R+ ) . Then 



Yim Lj{vn)F{vn,g) ^ lim / uj{x)F{x,g)dvri{x) ^ uj{v)F{v,g), 

n — >oo n — 'oc /q 

since x i— > lj{x)F{x, g) is continuous and bounded on the compact set containing the 
supports of the measures i/„. Since U!{v) is nonzero, it follows that F is continuous. 

D 

5.2. Invariance properties. 

Lemma 5.2. Suppose Xt converges on compacta to the quasistationary distribution 
(j), when started at the compactly supported distribution v . Let i„ be any sequence 
such that Ft^ has a continuous limit, as described in Lemma \5.1\ Then F^*"'{i', •) 
has density F'-*"\v,'E.'^)(j) with respect to Lebesgue measure. If convergence on com- 
pacta holds for a set of initial distributions which is dense in 9Jlf (M'*') (in particular, 
the distributions satisfying (IDC)), then F^*"''{i', •) has density _?"(*"' (jy, R+)(/) for all 
iyemi{R+). 

Proof. Convergence on compacta implies that for any z, z' > 0, 

Jo 4>{y)dy 
We may normalize cf) so that J„ (f>x{y)dy is either 1 or oo. Letting z' — > oo, we get 

(56) F(*")(i.,[0,z]) = i^(*")(i.,R+) r ,Px{y)dy 



JO 

for any v with a compactly supported density, and any appropriate sequence i„. 
By continuity, if we take a weakly convergent sequence of (IDC) distributions v, we 
see that (|56|l holds for all compactly supported i^. D 
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Lemma 5.3. Suppose Xt converges on compacta to the quasistationary distribu- 
tion (j), when started at the compactly supported distribution v. //_F'*"-'(z^, M+) = 1, 
then the limit lim„^oo '2i„('^, ?") exists, and is equal to L (/)(y)Py{ra > rjdy. More 
generally, if K :~ Ihiiz^oo k{z) exists, then the limit Ihiin^oo cit„{v,r) exists, re- 
gardless of the value o/_F'-*")(i/, R+) (recall the notation from (I2U|I ). and is equal to 

(57) a(*")(i/,r) = F^'-\i^,R+) f (t){y)Fy{T9 > r}dy + {l ~ F^'"^i',R+))e-^'- . 

Proof. By Lcnima l^^ pi^ri) (^jy^ .) j^^g density i^(*") (iy, R+)(/i^ with respect to Lebesgue 
measure. If g : M^ ^ R is any bounded continuous function such that 5(00) — 
hniy^oo g(y) exists, weak convergence imphes 

(58) lim / Ft^{v,dy)g{y)^ F'^'-\iy,dy)g{y) + {\ - F'^'-\v,M+))g{^). 

In general, 

«(*")(;., r) = I F^'"\y,dy)¥y{To > r} . 

If K has a Umit at 00, the function y t-^ Fy{Td > t} is continuous at 00 for each fixed 
r, from which H57|) follows. If F^*")(i/, R+) = 1, then the term at infinity vanishes 
in any case, so the convergence of k is irrelevant. D 

Proposition 5.4. Suppose the conditions (LP) and (GB), are satisfied. Let (i„) 
be any sequence such that Ft^ and tOt^ converge to continuous limits, as described 
in Lemma \5.1\ Then 

(i) The product uj^*"''{x)F^'^"'{x,M.'^) is X-invariant. It satisfies 

(59) c^(*")(x)F(*")(x,R+) = F(*")(l,R+)^4^. 

(ii) If K :— liniz^oo 't(^) exists, then oj^*^"' (x){l—F^'^"' {x,M.'^)) is K -invariant. 
If K is finite, then 

(60) Lo^'-\x){l - i^(*")(a;,M+)) - (l - ^(*")(1,R+))M4. 

■0k(1) 

If K ^00 then F(*")(a;,M+) = 1 for all x. 
(iii) If K := lim^^oo i^iz) exists and is finite, thenoj^'^"' is a convex combination 
of -0^ and ipK ■ Specifically, 

(61) a;(*")(.) = f (*")(1,M+)M^ + (1 - f(*")(i,k+))M^. 

(iv) If K := lim^^oo '^(2) exists, then for any compactly supported v and s > 0, 

(62) E^cj(*")(X,) =a(*")(iy,s)cj(*")(iy). 

(v) Either Xt escapes to infinity from every compactly supported initial distri- 
bution, or (f)\ is the density of a X-invariant distribution for Xt; that is, 

(63) P;0A = e-^VA- 
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Proof. Condition (GB) allows us to apply the Bounded Convergence Theorem, to 
show that for any positive s, 

Ei.w^ ^\X,) = E^ hm -— — — - 

(64) "^°° IriJTg > tn\ 



= ,„(*") 



LO 



{v) lim P^{rg >tn + s\TQ> i„}, 



which proves I)(i2|l . 
Define 

Then 



^xl{Xs) = E^ lim lim 



^— i-oo n — *oo 



= cj'^*"^(x) lim lim 
(65) 



iJTa > t„} 



2^oo n^oo P2:{T(3 > f„} 
/>oo 

= w(*")(x) lim lim / FtSx,dy)ry{X, < z] 

z—*oo n—*oo J r^ "■ 

/•OO 

= tj(*")(x) lim / F{x,dy)Ty{Xs<z\ 

/•OO 

= w(*")(x)i^(x,R+) / 0A(y)P4T9 > s}dy 

/>C30 

= 7(0::) / 0A(y)Py{ra > s}dy. 



/o 

Note that the step from the third to the fourth lines uses the fact that the function 
y I— > PyjXg < z} vanishes at cxd for each z, while the following step uses Lemma 

Suppose there is some initial distribution u satisfying (IDC), from which Xt does 
not escape to oo. Then Proposition 12.31 tells us that L (j)\{z)dz is finite. There is 
some sequence i„ such that F*^*"^(z^, M+) > 0. By the relation (|19|l . it follows that 
there must be some x such that i^((^^)(a;,K+) > 0. Given any positive s and t, 

E,7(X,+t)=E,Ejf,7(X,) 



E,7(X;) / ^x{y)^y{Ta > s}dy 
Jo 

/•CO /-oo 

7(.t) / cl>x{yWy{T9 > t} (t,x{y)Fy{Td>s}dy. 



Since this is also 

/•oo 

7(x) / c^x{y)Py{T9>s + t}dy, 
Jo 
it follows that the function 

/>oo 

(t)x{y)¥y{To > s}dy 

is multiplicative, and so may be written in the form e^''* for some rj. Observe now 
that 7(0;) > 0. Since 7 is an 77-invariant function, it must be in the domain of iL, 
and satisfy Lip = —ijip. By Lemma 12.21 and the fact that 7 is nonnegative and not 
identically 0, it must be that 77 < A. Lemma 17771 shows then that cpx is A-invariant, 
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proving (|63|1 . This also proves H59|l . in the case when Xt does not escape to infinity 
from every initial distribution. 

Suppose now that F(z^, R+) = for all (IDC) distributions v. By continuity, 
F{x,M.+) = for aU x. Then ^ holds trivially. 

It remains only to show (|60() — equation H61() merely combines this with (|59|l . 
We repeat the computation of 1)65(1 . now with 

^(x) := Lu^'"\x){l ~ F(*")(x,R+)). 

We have 



z — >-cx3 n — ^■cx: 



(66) 



E^7(X^) = E^ lim hm 



l{Td > tn} 

lim / F<^'-\x,dy)Py{Xs>z} 



F(*")(x,M+)) lim P„|Xs >z| 



^{x)e' 



-Ks 



In the last step, we have used the fact that oo is a natural boundary, which implies 
that 

lim Vy{Xs > z| = lim P„{ra > s]. 

a 

Corollary 5.5. Suppose (LP) and (GB) hold. Let tn -^ oo he a sequence such 
that Ft^ and ojt^ have continuous limits, as described in Lemma \5.1l Then either 
i^'*"^(i/, R''") — for all compactly supported v or F^*"'{iy', M"*") > for all compactly 
supported v. If K — lim^^oo 't(^) exists, then either i^''*"'(z^, M'^) = 1 for all 
compactly supported v or i^^*"'(i/, R^) < 1 for all compactly supported v. 

Proof. Suppose there is some ly such that F(*"'(z^, R+) = 0. By H19() . there must be 
some x such that F(*"'(a;,R+) = 0. By ProDOsition l5.4l we know that 

E,c^(*")(X,)i^(*")(X„M+) = 0. 

By regularity, and the fact that uj is positive, it follows that F'^*"'(a;,R+) — for 
Lebesgue-almost-every x. Suppose there were exceptional x' such that F^*'"\x' , M+) > 
0. Then 

E,,w(*")(X,)f(*")(X„R+)>0, 
which is impossible, since the distribution of Xg has a density with respect to 
Lebesgue measure. From (|19() . it follows that F*^*")(z/',R+) = for all compactly 
supported v' . An identical argument, using ((60(1 in place of 1(59(1 . shows the second 
part. D 

Proposition 5.6. Suppose the conditions (LP) and (GB) hold, and that K := 
lim^^oo k(-z) exists, and is not equal to A. Then the following are equivalent: 

(a) For any sequence f„ such that Ft^ and uot^ converge, _F'^*"-'(i^, R+) is con- 
stant in V . 

(b) If V is any compactly supported initial distribution, then limj^oo -Fi('^, R'*') 
exists. For each v, the limit is either or 1. 

(c) There is some v e ^^(-l^^) smc/i that linif^oo -Ft(i^, R^) exists, and is 
either or 1. 
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(d) For all x G R+, the limit linit^oo wt(a;, M"'') exists. For each x, the limit is 
either ipK or '0a- 

Proof, (a)^(b) Suppose that F*^*"^(;y, R+) is constant. By Proposition 15.41 we 
know that Ci;^*"-'(x) is either A-invariant or iiT- invariant, since F(*"^(a;,R"'") and 
1 — F(*")(a;,K"'") cannot both be 0. But one of these must be 0, since ti;'^*"^(a;) 
cannot be both A-invariant and i^-invariant. Thus F'*"-'(i/, R"*") is either or 1. 

Suppose, now, that there were a sequence i„ such that F'^*"'(z^, R+) = 0, and 
another sequence t'^ such that F(*")(z^,R+) = 1. Since Ft{x,-) is continuous in i, 
it would be possible to find an intermediate sequence t" such that the Hmits i^(*") 
and w^*"^ exist, and i^(*"^(a;,R+) G (0, f). This proves that i^(a;,R+) is sequence- 
independent. 



(b)^(a) From the invariance properties of Proposition lS^ it follows that F{x, 
is constant in x. Suppose now that limt__+oo ^t(2;,R^) = 0, and let i^ be any com- 
pactly supported probability. Choose e > and z > 0. For each x, we may find 
t{x) such that Ft{x, [0, z]) < e for all t > t{x). Thus 

r I -, \'^x{Xt < z}di^(x) r , . -, 

MX^<^\ro>t}-^j^^^^^f^<H^:tix)>t}+e, 

which converges to e as i —> oo. Letting e ^ proves that Ft{i^, [0, z]) = 0. A 
similar proof holds if F{x, M+) = 1. 

(b)=^(c) This is trivial. 

(c)^(a) Suppose first that limt^^ Ft{i^,R+) = 0. Then F(*")(i^,R+) = for 
any sequence i„. By Corollarv l5.5l it follows that F*^*")(i/',R+) = for any v' . An 
equivalent proof works if limj^oo Ft{v, R^) = 1. 

(b)<^(d) This is a direct consequence of JHTJ. D 

5.3. The Dichotomy: Proving Theorem 13. 31 for A < liminf.^^^ k(z). Assume 
A < liminf^^too k{z). If L (f)x{y)dy = oo, then the process escapes to infinity from 
every starting distribution satisfying (IDC), by Proposition l2.3l Since (IDC) distri- 
butions are dense in 9Jlf (R+), by the continuity of F we see that i^(*"^(z^, R+) = 
for any compactly supported v. 

Assume now that J^ (f)\{y)dy is finite. We need to show that 

(67) hm hm g^L^ = 

(We drop the compactly supported starting measure from the notation.) 

We first prove H67() with t running only over integers. For positive integers n and 
positive real z, we define 

(68) /(.,n):^^P">t 

Choose Zq such that K := infj^>2j, k(t/) > A. Observe that for z > zq 
P{X„+i > z} < P{X„ < z and X„+i > z} 

(69) + P{X,i > z and X„+i > z and Xt > zq Vt G [n, n + 1]} 

+ F{Xn > z and Xn+i > z and 3i G [n,n+ 1] s.t. Xt < zq}. 
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Thus 

f{z,n+l)<F{Xn+i> z\Xr^<z} 

(70) +/(^,n) J^\;^"~f\ [P{X„+1 >zkXt>zo,yte [n,n+l]\Xn>z} 

+ P{X„+i > z&3t e [n,n+l] s.t. X* < zq | X„ > z} . 

By the strong Markov property (taking the first time the process hits the level zo 
after time n), the last term in brackets is bounded above by 

(71) P,„{3ie[0,l]s.t. Xt>z}. 

The middle term in the brackets in (|70|l is bounded by e~^ , since the process is 
being killed at rate at least K while it is above zq. 
Substituting these bounds into (|70() . we get 

/(z,n+l) 

V{Xn<z} 



(72) 



"{^n+l < z} 



-K 



{X„+i > z\Xn<z] 

,„{3ie [0,1] s.t. Xt>z])f{z,n) 
We write this in the form 

/(z, 71 + 1) < c(z, n)J{z, n) + 7(2, n), 
where 

(73) c(z,,i):-(P,„{30<i<l : X, > z} + g-^') (^^^^i^^ly), 

and 

(74) 



7(z, n) := P{X„+i > z I X„ < z} (- 



"{X„<Z} N-l 



'{^«+i < z] 



From Lemma 17.51 we know that 



lim 



'{Xn < z} 



n^^ P{X„+i < z} 

Substituting into H73|l yields, 

limsuplimsupc(z,n) < e*^-^-^""'^ + e- lim sup lim sup / Pzo{30 < t < 1 s.t. Xt > z] 

z — >oo n — >oo z — ^00 n — >oc> JQ 

(The second term on the right is because 00 is an inaccessible boundary.) 
By quasistationary convergence on compacta, 

lim sup lim sup 7(2,71) < lim sup / Pj,{Xi > z^d^z{y)i 

z — >oo n — ^00 z — >oo Jo 

where /i^ is the distribution on [0,z] with density proportional to (jjx- Since 00 
is inaccessible, it follows that lim^^oo Pyj-'^i > z} = for each fixed y, so that 
limsup^^oo l™sup„_^oQ 7(2, n) = 0. Lemma 17.61 implies then that 

lim sup lim sup /(z, n) = 0. 



z — ^cx3 n — >C30 
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Now we need to extend the proof of l|H7| to real times t. We have 

y y H^t > z] 

hm hm sup — ? f 

(75) I "^ - J 

, , F{3te[n,n+l)s.t. Xt> z] P{X„ < 2} 

< hm hm sup —7 r sup —7 r- . 

- z^oo n^J^ P{X„ < Z} tG[«,n+l) P{Xt < Z} 

The numerator of the first term on the right is bounded above by 

F{Xn > z}+F{Xn < z} -Plate [n,n+l) s.t. Xt> z\Xn<z}. 

By quasistationary convergence on compacta, the hmit as n ^- 00 of the right-hand 
term is 

<j>x{x)dxj I Fy{3t e [0, 1) s.t. Xt > z]c^x{y)dy. 

Since cx) is inaccessible, this goes to as z ^ 00. 
The expression in H75|) is then bounded by 

p{x„<4P{x„>z} 

lim sup lim sup sup ^, ,^ ^ -, t- , 

,^J^ n^J^ t^[n,n+l] H^t < z} P{X„ < z} 

which we have already shown to be 0. 

5.4. The Dichotomy: Proving Theorem l3.3l for A > limsup.. _, k(z). Heuris- 
tically, the proof is very simple, though the notation may be confusing. Consider 
some large z, large enough that the killing rate is always smaller than A above z; 
and choose some z' > z, large enough that the probability of moving from above z' 
to below z in one unit of time is small. Eventually, mass below z must be declining 
at rate A. The mass above z' is being killed at a strictly lower rate. It could be, 
then, that the mass above z' blows up relative to the mass down below, in which 
case we have escape to infinity. Otherwise, the mass above z' must be sliding down. 
Can this happen? The obstruction is the narrow pipe between z and z'. In one 
unit of time, only a small fraction of the mass above z' can make its way down 
below z. Thus, if mass is escaping downward, it must be spending some time in 
[z, z'\. Convergence to a probability distribution on compacta implies, though, that 
eventually the mass on [z, z'\ — taking up a significant fraction of mass from [z', 00) 
— will be a small fraction of the mass on [0, z]. Putting all of this together, we 
conclude that the mass above z' must be much smaller than that below z. 

We now translate the heuristic into formal mathematics. Fix z large enough that 

supk(2/) < {K^X)I2=:K, 

y>z 

and define 

._ XT (t^>S~y)'^y 

Jo (t>\iy)dy ' 
Associate to z some z'(z) > z (assumed increasing in z) such that 

P^'{3te [0,1] : Xt<z} <e. 
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(This is possible, since we have assumed that oo is a natural boundary.) We define 

a„:^P{X„<z}, 
6„:=P{z<X„ <z'}, 
c„:=P{X„>z'}. 
We have the following relations: 

c„+i>e-^c„(l-P,.{3ie [0,1] : X* < z}) - 6„+i, 

bn 

hmsup — < e, 

n — ^oo O^n 

lim = e -. 



Thus, 



>- e (1-e) — 



Set i? = e^"^(l - e)2. We may find N such that for all n> N, 

>R^- 2e. 



Suppose z large enough so that R> 1. Then cither 
iim — — CO or lim sup — < 



These translate into 

2e 



lim f{z,n)—oo or limsup/(z, n) < 



i?-l 

for z sufficiently large, where f(z,n) is as defined in (|68|l . Note that f{z,n) is 
monotone decreasing as a function of z, so that if the second limit holds for any z, 
it holds for all higher z. Sending z to oo makes e — > 0, so that cither 

(76) lim f{z,n) — oo for all z or lim lim sup /(z,n) = 0. 

n — >oo z — >oo ^ ^QQ 

The same argument as in section 1)5.31 shows that the second condition in H76|) is 
equivalent to 

lim limsup/(z,i) = 0. 

Similarly, if lim„^oo f{z, n) = oo for all z, we have for any z' > z, 

,. . rr, ,,^r . ,P{X„>z'}-P{X„>z'&3te (n,n + l) : Xt<z} 
limmf f (z,i) > hmmf — ^ - 

t^oo -"^ ' ' - n^oo P{X„ < z'} + ¥{Xn > z' & 3t £ (n, « + 1) : Xt< z} 

^^. . f(z',n)(l-P{3ie (n,n + l) : Xt < z I X„ > z'}) 

> lim mt ^ — -— 

- n-.oo 1 + f[z', n)P{3i e (n, n + 1) : X* < z | X„ > z'} 

= liminfP{3te {n,n+l) : Xt<z\X,,> z'y^ - 1. 

Since oo is a natural boundary, this diverges to oo as z' — > oo. Thus, we have shown 
that 

(77) lim /(z,i) = oo for all z or lim limsup/(z,f) — 0. 
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If we were assuming (GB), we could apply Proposition 15.61 here. Instead, we 
note that v i— > i^(*"'(z^, R+) is continuous. Since we have shown that it takes on 
only the values and 1, it must be constant. The case F(*"'(z^, IR+) = is the 
same as escape to infinity; the case F'^*")(i/,M+) = 1 is (by Lemma f5.2|) the same 
as convergence to the quasistationary distribution 0a. 

5.5. Proof of Theorem l3.4L If A = K, then Lemma lsm together with (|6fl) implies 
that every sequence t^ ^ oo has a subsequence i„ such that lim„_^oo ij-'t„ (x) = 
il}\{x). Thus, ip\ is the limit. The same argument applies to at„{x,r), proving the 
claim about the killing rate. 

Suppose now that X < K. Since A is the bottom of the spectrum of iL, it 
follows (as we have made explicit in Lemma I2.2|l that there is some a; > for 
which tpK{x) = 0. Since w is strictly positive (except perhaps at 0), it follows 
from H6UI) that F(*")(x,E'^) = 1 for any sequence (t„) for which the limits exist 
and are continuous. By Proposition l5.6l limt^oo Ft{i',M.'^) — 1 for every compactly 
supported v. The results of Theorem 13.41 follow then from Lemma 15.21 and Lemma 

E31 

Suppose, finally, that X > K. If Xt escapes to infinity from any compactly 
supported initial distribution, it follows, by Proposition 15.61 that Xt escapes to 
infinity from every compactly supported initial distribution. Suppose, then, that Xt 
does not escape to infinity from any compactly supported initial distribution. Then 
/o 4'\{y)dy < c», so we may apply Theorem 13. 31 to see that limj^oo i^((j/, M+) = 1 
if v satisfies (IDC). By Proposition l5.6l this holds then for all compactly supported 
v. The result then follows from Lemma [5. 21 

6. An example 

One Markov mortality model that has received considerable attention is a discrete- 
state-space process that was first defined by H. Le Bras |Le 76j . (A somewhat sim- 
ilar process appeared in the mathematical literature, presumably independently, in 
|Cav78| .') In this process the states are positive integers, and the motion is only 
increasing by single steps. The rate of transition to the next higher state, and the 
mortality rate, are both directly proportional to the current state. 

We define a continuous-state-space analogue to this process on [l,oo) by the 
stochastic differential equation 

(78) dXt = aXtdWt + bXtdt, 

where cr is a positive constant and b is a constant larger than a^ /2. The process 
starts at Xq = 1, is killed at the rate kXt, and is refiected when it hits 1. By Ito's 
formula we see that Xt has the same dynamics in (1, oo) as the geometric Brownian 
motion 

At=exp{a(I¥t+k)}, 

where h — ^ — §• Equivalcntly, then, we could consider the Brownian motion with 
drift: 

Yt^Wt^ It, 

killed at a rate ke'^^ and refiected at 0. This will have the same mortality distri- 
bution as Xt- We have B{x) = bx, where b is positive. The left boundary is 
regular, while the right boundary (X) is natural. Since k{x) ^ oo as x — ^ (X), we 
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can apply Theorem 13.31 This guarantees the convergence to a quasistationary dis- 
tribution from any initial distribution with a compactly supported density, as long 
as /„ (j)\{y)dy < oo. Condition (GB) is satisfied, as we discuss in the comments 
preceding Lemma 12.51 since the drift is constant and k is increasing. This allows 
us to apply Theorem 13.41 to conclude that L (t)\{y)dy < oo. (This would follow 
from our computation in any case.) 

Eigenfunctions 4>x for the generator of Yt are given (up to a scaling constant) by 

where ^\ : [1, cx)) -^ M satisfies 

(79) ia2(a;2eA(x))" - bixU^))' - kxU^) = -XU^), 
with the boundary condition 

(80) eUl) = (^ - ^^(1)- 

In a sense, ^\ is the function of interest, since this is the density of the quasista- 
tionary distribution in the original coordinates. Formally, though, we need to refer 
to (j)x, since our theorems are stated in this normalization. In particular, it is (fix 
that needs to be integrable.) 

If we define a new function ip by 

we get 

zV(^) + zi^'iz) - {v^ + z^)i'{z) = 0, 
where 

(81) 

(Note: In order to find the transformation, we begin by writing S,x{x) = x''V'(\/aa;), 
and solve for a and r.) This is the equation satisfied by the modified Bessel func- 
tions. Thus the general solution for (j) is 

where v is given as a function of A by H81|l . On any compact set of z and A, S,x{z) is 
uniformly continuous in z and A (by Theorem 1.7.5 of jCL55| ). which implies that 
Ci and C2 are continuous. 

We need to compute A. By Lemma 12.21 we know that ^x is positive, but ^x 
changes sign for all A > A As 2: ^ 00, Iu{z) is increasing like e^/^/z and Ki, 
decreasing like e~^ j ^fz. Togther with the continuity of the coefficients, this implies 
that Ci(A) = 0. The top eigenfunction that we seek is a multiple of ly. 

The transformation of t,x to V' shifts the boundary from 1 to 






xq = y^k/a. 
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We have 

f (a;) / b 



(82) 



<f>{x) 



t2 




Thus, the boundary condition turns into Kl{xo) = 0. 

There cannot be a solution to (|82|l for i/ real and positive, since then K^{xo) is 
positive, while the derivative is negative. By the integral representation (9.6.24) of 

/•OO 

Jo 

/>OC 

Kly{x) = - / 6-^^°*^'^* cosh(i) cos{yt)dt. 

Observe, first, that Kiy{x) has infinitely many zeros close to a; = 0, but also a final 
zero; the location of this final zero is monotonically increasing in y, and unbounded. 
The same is true for the derivative K[ [x) , and the final zero of K' is larger than 
the final zero of K. The zeros of K and K' are interleaved. For jy = 0, Kiy{xo) is 
positive, while K-y{xo) is negative. For y close to this ratio remains negative. As 
y increases, eventually the final zero of K'^ reaches xq, at j/ = y{xo)- (Remember 
that a;o is itself a function of the process parameters a and k.) The limiting rate of 
killing is then 



a^ 



A-A(iy) = — 



(^-1) ^e 




The quasistationary distribution is a constant multiple of 

This function behaves asymptotically as 

as x ^ oo, according to formula (9.7.2) of |AS65j . We see that (t>\_{x) = <7e'^^^\{e'^^) 
is integrable, confirming the conclusion of Theorem l3.4l 

7. Technical lemmas 
7.1. Facts about generators and resolvents. 

Lemma 7.1. The adjoint semigroup for the diffusion with killing is generated by 
the operator with domain 'D = {/ G L^ n C : C* f G L^}, acting as L* . 

Proof. The Feynman-Kac construction already defines a strongly continuous con- 
traction semigroup and an adjoint strongly continuous contraction semigroup. Ini- 
tially, the adjoint semigroup acts on 'BV([0, r]). We observe, first, that the image of 
any measure under the resolvent has no atoms at the endpoints, since we have as- 
sumed instantaneous reflection (in the case of regular boundaries) or instantaneous 
entrance (in the case of entrance boundaries). Now, let i?^ be the resolvent for the 
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semigroup with internal killing, and i?^* the adjoint resolvent with internal killing. 
For V any finite measure on IR+ and A C M+, 

/•oo 

i?XA) = E^ / e""*l{X( e A}dt, and 
Jo 

/•oo 

We see then that R*^v and B^*v are equivalent measures. Since we know from 
Theorem 13.3 of |Fel52| that R*^v has a density, so must B^*v have a density. Both 
are bounded in total variation by 1/a. This tells us, in particular, that all of L}- is 
mapped into i^ by the resolvent. Since the range of the resolvent is dense in the 
range of the semigroup, this means that the restriction of the adjoint semigroup 
action to V- is a strongly continuous contraction semigroup on L^. 

Call the generator of this adjoint semigroup with internal killing fJ*. Formally, 
the adjoint generator must act as the formal adjoint of the generator, namely, as 
£*, on functions in the domain which are twice differentiable. It remains to show 
that 2) is a core for fJ*. By Proposition 1.3.1 of EK8^, it suffices to show that 2) 
is a dense subset of the range of the resolvent, and that the resolvent maps D into 
itself. 

We first need to show that 1) is actually contained in the range of the resolvent. 
That is, for every f G T) we must find g G L^ such that R'^*g ~ f ■ Equation 6 
of section 4.11 of |IM65j represents the a-Green's function (which is the resolvent 
density) for the killed diffusion by 

(83) Sa(y,a;) = g„(x,2/) = <^ , s b{x^ ■( ^ 

where gi and (72 are increasing and decreasing eigenfunctions, respectively, for L, 
and c is a constant. Define g — af — L* f and h = R^g. By the definition of T), f 
is in L^ and twice continuously differentiable. Thus g is continuous, and it follows, 
by Lemma [7.21 that h is continuously twice differentiable, so that fi* acts as the 
differential operator L* on h. This implies that both / and h are solutions to the 
differential equation 

au — L*u = g, 
so that 

Ii*[f-h]=a[f~h]. 

We show first that h satisfies the boundary conditions 0) at 0. The boundary 
condition for the image of i?^ is given on page 131 of IIM65I as 

(84) {1-po) lim>(a;) = ipo lim e^ (-),/,' (x). 
This is equivalent to the equation for the Green's function 

(85) Vy, (l-po)limG„(a:,y) = ipolime^(-)^%^. 

x^o 2 x^o ax 

Now, from the representation of the Green's function as H83|) . we see that Ga{x, y) = 
eB(v)-B{x)Q*^(^^^yy This means that 

-ldG*Jx,y) 



lime^(-)^^4^=e^(^'lim 

a;— >0 OX x~^0 



2 dx 



h{x)Ga{x,y) 
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This translates into the condition on h = R'^*g of 

(1 - po) hm e-^(^)/i'(a;) = po 1™ \h\x) - 2b{x)h\x)] , 

which is exactly the condition in Q. Similarly, if r is regular, h satisfies the 
condition |(SJ), and if r is an entrance boundary h satisfies ©• Thus f — h is a.n L^ 
solution to £j*u = au, satisfying these boundary conditions. We show in Lemma 
17^ (copying |Fel52| ') that this implies that f -h^O. 

We now show that i?^* maps T) to itself. We have just shown that the adjoint 
resolvent takes all continuous L^ functions to continuously twice-differentiable func- 
tions which satisfy the appropriate boundary conditions. If / = Ra*9i where g G T), 
then f E G, and 

L*f^af-geL\ 
Finally, if / is in the domain of fi*, we can represent / as i?^*5 for some g E L^, 
then we can represent g as a limit of continuous functions gi. Then fi — R^gi is 
in D, and satisfies \\fi — f\\ < a^^Wgi — g\\, implying that limi^oo ft = f ■ As we 
have already shown that the functions gi are in B, this shows that T) is dense in 
the domain, completing the proof. D 

Lemma 7.2. Letting G*^{x,y) he the adjoint Green's function given in (|83|l . and g 
any continuous L^ function, then 



h{x) := / Gl{x,y)g{y)dy 
Jo 

is continuously twice differentiable. 

Proof. The technical center of this result is a merely tedious calculus fact, based on 
the differentiability of G* away from the diagonal. We isolate the details as Lemma 
17. HI It remains only to show that the conditions of that lemma are satisfied by G* 
for n = 2. 

Since we assumed that b is continuously differentiable, and that k is continuous, 
it follows that B is twice continuously differentiable, as are the eigenfunctions gi 
and 32 of A*. We have G*{x, y) = gi{x A y)g2{x V y)e^(^^. We have 

''^^'^'^^ \c{g,{x)e^(^y)"gM iiy>x. 

Since gi is increasing and g2 decreasing, the factor which depends on y is bounded. 
This implies that G2{x, y) and Gl{x, y) are bounded in y, and that their continuity 
and convergence are uniform in y. As for G^ , we have 

Gl'^ix) = (5i(x)e^(-))'52(2/) - {g2{x)e^^''^)' g^{y), 
which is continuous. Since G* is continuous, it follows that Gg is identically 0. D 

Lemma 7.3. LetT(x,y) be a real-valued function on [0,r] x [0,r], which is n-times 
differentiable as a function of x for every interior point y ^ x, where n is a positive 
integer. Let g € L^(0,r) be continuous, and define 

h{x) := / T{x,y)g{y)dy. 
Jo 

Define Tk{x,y) — d^T /dx^{x,y) for k <n, and let 

^kiy) :=limrfe(a;,y) - limrfe(a;, y). 

xiy x\y 
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(87) 



Suppose that 

• Tk{x,y) is bounded as a function ofy, for k < n. That is, 
(86) sup\rk{x,y)\ < oo 

• Tk{x, y) is continuous in x, uniformly in y ^ x, for < k < n. That is, 
lim sup \rk{x + e,y) — erk{x,y)\ — 0, and 

e^0+ ye{a,x)u{x+e,r) 

lim sup \Tkix- e,y) ~ eTkix,y)\=0. 

e— >0+ j/g(0,a;-e)u(a;,r) 

• The convergence to the k-th derivative is uniform in y ^ x, for I < k < n. 
That is, 

lim sup e^^\rk-i{x + e,y) -Tk-iix,y) - eTn{x,y)\ ^0, and 

6^0+ ye(a^x)u{x+e,r) 

lim sup e^'^\Tk-i{x -e,y) -Tk-i{x,y) -eTkix,y)\^0. 

£-^0+ ye(Q,x-e)u{x,r) 

• r„_j^ is continuous, and V^ is zero for fc < n — 2. 

Then h is n-times continuously differ entiable on (0,r). The n-th derivative is 
given by 



(89) 



hnix):=fix)T^_,ix)+ / T„ix,y)fiy)dy. 
Jo 



Proof The proof is by induction on n. First take n = 1. We have assumed that 
Fq is continuous, and Fi is continuous in x, uniformly in y. Thus hi exists, and is 
continuous. We need to show that for each x E (0, r), 

= lim e^^ [h{x + e) - h{x) - ehi{x)] 



= lim e 



9(y){^(x + e,y)-T{x,y) -eTi{x,y)}dy - eg{x)r^{x) 



We consider limit from the right; the limit from the left is essentially the same. For 
y outside of [x, x + e], the quantity in braces is o(e), uniformly in y. This leaves us 
with the task of showing that 



= lime'i 
We can rewrite 



giy) {^(x + e,y) -T{x,y) ~ eTi{x,y)} dy - eg{x)TQ{x) 



rV rx+e 

T{x + e,y) -r{x,y)^ I Fi(z,y)rfz + F^(j/) + / Ti{z,y)dz 



a;+e 



Ti{z,y)dz + T^{y), 



transforming the limit into 



a^+e 



9{y) 



= lime"i 

40 

It will suffice, then, to show that 



x-\-e 



{Fi(z,2/)-Fi(x,2/)}dz + r^{y) 



- 9ix)Toix) dy 



= lim sup g{y) 



40 



yG[x,x+e] 



{ri(z,y) -Fi(a:;,y)}(iz + T^iy) 



gix)T^ix). 
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By the assumption (|86|l . boundedness of Fi, the integral is 0(e). Since g and T^ 
are continuous, the remaining difference goes to with e. This completes the proof 
for n = 1. 

Now, take n > 2, assume the lemma true for n — 1, and that the conditions 
are satisfied for n. If we substitute Fi for F, the conditions are satisfied for n — 1. 
Furthermore, since Fg = 0, we have 



h'{x)^hi{x)= Ti{x,y)g{y)dy. 
Jo 

It follows from the induction hypothesis that h'(x) is {n — l)-times continuously 

differentiable, which completes the induction. D 

The following lemma is copied from the results of sections 9 through 11 of |Fel52| . 
adapted for the case of nonzero internal killing. With only a few small exceptions, 
the arguments are identical, though we have specialized them somewhat for the 
problem at hand. 

Lemma 7.4. For a > 0, if 4> E L^(0,r) is a solution to the differential equation 
L*(j) — a(j) satisfying the boundary conditions Q, and ^ if r is regular, then 

Proof. Let u be any solution to the forward eigenvalue problem Lu = au on (0, r). 
We note first that u cannot have any positive interior maxima. If a; is a local max- 
imum, then > u"{x) = 2(K(a:) + a)u{x), which imphes that u{x) < 0. Similarly, 
there is no negative minimum. It follows that u has limits at both endpoints. It 
also follows that two solutions ui and U2 with ui{xo) = U2{xo) = 1 do not meet 
at any other point. (Recall that xq is an arbitrary point in (0,r).) Thus, for each 
s £ (0,r) \ {xq}, there is exactly one solution u such that u{s) = and u{xo) = 1. 
Letting s — > r, the derivative u'(0) converges to a negative limit, yielding a solution 
which is positive and decreasing on {xo,r). 
If V is any solution to Lv — au, 

(90) w'(x) =e--^(^)|w'(a;o) + 2 / {a + K{z))v{z)dz 

y J Xq j 

A consequence is that v'{x)e^^^^ always has a limit at both boundaries. 

Let us take a solution with v(xq) — and v' {xq) > 0. All terms in (|90|l are 
positive for x > xq, so v is nondecreasing. It follows that 

< v'{x) < w'(xo)e--^(=") + 2(a + K)v{x)e-^^^'^ f e^^'Uz =: Bi{x) + v{x)B2{x). 

J xa 

Suppose r is regular. Since and r are regular, Bi and B2 are integrable on 
(0, r). It follows that v is bounded on (0, r). From this w, and the nonincreasing u 
mentioned earlier, we can compose two positive nonincreasing solutions mq and ui 
on all of (0, r), such that 

limuo(a;)=0. Van e^^'-'^'^ix) = -I, 

X — yr X — >r 

lim ^1(2:) = 1 and lim e^^^^u'i{x) = 0. 



Both solutions will be positive at 0. As has already been pointed out, the functions 
Vi :— e^Ui are eigenfunctions for the adjoint operator, satisfying L*Vi = avi. Since 
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they are linearly independent, they span the solution space, and we can represent 
4> as aofo + aivi for some real numbers ag and ai. By (jSJ, if ao and ai are not both 
zero, 

> lim ^"'^^^^(^) -^ 



i'/''! 



x) — b{x)(j){x) — flQ 



This implies that ai and ao have the same sign, or ai — 0. Without loss of 
generality, suppose both coefhcients are positive. We then have 

< lim e"''^"V('^) ^ aoUoiO) + aiuijO) 

- x^o ^(jj'ix) - b{x)(j){x) aoe^(^)u[|(0) + aie-S(^)u;(0) ' 

Here the numerator is positive, and the denominator negative, which is a contra- 
diction, proving that ao and ai are indeed both zero. 

Consider now the situation when r = cx3 is a natural boundary. Let u be an 
eigenfunction as above, with u(xo) > and u'{xo) > 0. By (|9()|l . 

u'{x) > ae-«(^) / e^'-^^dz. 

J Xq 

As oo is a natural boundary, it must be that e^'^'-^' /^ e^^^-'dz is not integrable. 
Together with the fact that u is increasing, this implies that linij^^oo "(a;) = cxd. 
Furthermore, for x> xi > xq we have /3 > such that 

J X\ 

Again, since oo is a natural boundary, it follows that 

e^'''^^u{x)dx = oo. 



Since g L"'^, it must be that u = e^ is positive and nonincreasing, negative and 
nondecreasing, or identically 0. By the same computation as above, either of the 
first two cases implies that 

lim —r-. \ , , < 0, 

x^o (t>'{x) - 2h{x)(t){x) - 

contradicting the assumed boundary condition |0J. Thus = 0. 

Finally, suppose r is an entrance boundary. By the same argument as above, we 
see that solutions u passing through are unbounded at r, and when u'{xq) > 
for u{xi^) — 0, we have lim^j^r e^^^-'M'(a;) > 0. It follows that, for (j) to satisfy the 
boundary condition © it must be a multiple of e^uo, where u is a nonincreasing 
solution. As before, though, a nonincreasing solution cannot satisfy the boundary 
condition Q at 0. D 

7.2. Facts used in the proof of Theorem 13.31 

Lemma 7.5. Suppose (LP) holds. Then for any distribution v satisfying (IDC), 
and any positive z and s, 

(91) lim ^;^f;^- ^ f = e-A^ 

^ ^ t^^ V,{Xt<z} 

The limit is uniform in s. 



QUASISTATIONARY DISTRIBUTIONS 



37 



Proof. Using the notation of section WaI we need to show that 



(92) 

uniformly in s. Let 



lim 



-As 



t^oo x(t,5^[o,2]!; 



h{\) = sup(e-^'' - e-^") < T- - 1- 
s A 



For all t large enough that I(i,5C[o.2]) > Oi 



sup 

sG[0,oo) 



X(t + S,g^[0,; 



^it^9i[o,z]) 



-e-^^ 



sup 

se[o,oo 



1(^ + 5,5^(0,.]) At,9S,[o^z]e-n 



< 






1{t,9S,[o,z]) 



2^(^,5^(0,2]) 
An application of (|49|l proves that this converges to as t ^ oo. D 

Lemma 7.6. Let f.^c^-f : M+ x Z+ -^ M+ be maps satisfying 

(93) f{z, 71 + 1) < c(z, n)f{z, n) + 7(2, n), 

and c := lim sup^^^^ lim sup^^^^ c(z, n) < 1 . Lei 7 := lim sup^^^^ lim sup^^^^ 7(2, n) 
Then 

7 



(94) 



lim sup lim sup /(z,n) < 



z — ^cxD n — ^00 (^t Cj 

Proof. Choose e e (0, 1 — c), and find no and zq such that for any z > z^ and 
n > no, c(z,n) < c + e and 7(2, n) < 7 + e. Then for any z > z^ and n > no 
repeated application of l|93|) yields 

/(^,n)<(c + e)"-"°/(^,^o) + (7 + e)(l + (c + e) + --- + (c + e)"-"''-i) 

1- (c + e)"-"o 



(c + e)"-"«/(z,n) + (7 + 6)- 



1-c-e 



Thus 



lim sup lim sup /(z, n) < 



z — >oo n — *oo 



7 + e 
1-c-e' 



Since e was arbitrary, this completes the proof. 



D 



Lemma 7.7. ///q 4i\{y)dy is finite, then (j)\ is X-subinvariant for the semigroup 
P* . That is, 

(95) Pr(</'A)<0AP;^-l<e-^«</.. 

(96) / Py{To > s]<t>x{y)dy > e 

then 4)\ is X-invariant. 

Proof. Define the h-transformed operator 

.-..;.i-c,/*.)4g.(|-.)|-i/. 



As 
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This defines a diffusion on M+ with constant kiUing at rate A. If we let P, - be 
the corresponding semigroup, we have 

(This is an inequahty rather than equahty because the diffusion could explode.) If 
(|9(j|l holds as well, then 

/>oo 

[P:My) - e^-'0(y)] dy>0> P:^x{v)dy - e'^'. 

'0 Jo 

Since the right-hand side and left-hand side are equal, it must be that the integrand 
is zero almost everywhere. We conclude by continuity that P*4>\{y) = e~-''(j){y) for 
all y. D 

Note that this is essentially a special case of Theorem 4.8.5 of |Pin95| . 
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